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ABSTRACT 

We present results from the Chandra X-ray observation of Abell 13, a galaxy cluster that contains an 
unusual noncentral radio source, also known as a radio relic. This is the first pointed X-ray observation 
of Abell 13, providing a more sensitive study of the properties of the X-ray gas. The X-ray emission 
from Abell 13 is extended to the northwest of the X-ray peak and shows substructure in dicative of 
a rece nt merger event. The cluster X-ray emission is centered on the bright galaxy H of ISlee et al.l 
(|2001h . We find no evidence for a cooling flow in the cluster. A knot of excess X-ray emission is 
coincident with the other bright elliptical galaxy F. This knot of emission has properties similar to 
the enhanced emission associated with the large galaxies in the Coma cluster. 

With these Chandra data we are able to compare the properties of the hot X-ray gas with those 
of the radio relic from VLA data, to study the interaction of the X-ray gas with the radio emitting 
electrons. Our results suggest that the radio relic is associated with cooler gas in the cluster. We 
suggest two explanations for the coincidence of the cooler gas and radio source. First, the gas may 
have been uplifted by the radio relic from the cluster core. Alternatively, the relic and cool gas may 
have been displaced from the central galaxy during the cluster merger event. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (A13) — X-rays: galaxies: 
clusters — radio continuum: galaxies — intergalactic medium 



1. INTRODUCTION 

The high spatial resolution of the Chandra X-ray Ob- 
servatory allows us to study the intracluster medium 
(ICM) in clusters of galaxies in unprecedented detail. 
One interesting result is the interplay between the X- 
ray emitting gas and non-thermal radio emission. Ob- 
servations of several galaxy clusters have found X- 
ray cavities in the ICM coincid ent with bright radio 
emission (e.g., Pe rseus (A 426); iBohringer etai] 119931 : 
iFabian et al.l l2000h . The central radio sources produce 
the bubbles of non-thermal plasma responsible for the 
radio emission. These bubbles then rise, possibly by 
buoyancy, displacing the X-ray emitting gas and caus- 
ing cavities in the latter. 

A number of galaxy clusters contain large-scale diffuse 
radio sources not associated with active galactic nuclei. 
These sources are classified into two groups: radio halos 
which are centrally located, and radio relics. Radio relics 
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are extended diffuse radio sources often located at the 
periph ery of clusters of galaxies. iGiovannini fc Ferettil 
(2004) present a classification scheme for radio relics 
based on their obse rvational properties (but see also 
iKempner et al.ll2~004h . Radio relics typically have steep 
radio spectra ( a > 1 assuming S u oc y~ a \ and are lin- 
early polarized (|Giovannini fc Ferettil 12004 ). The X-ray 
surface brightness around most radio relics is low, making 
it difficult to study the relationship between the X-ray 
gas and radio emission. 

Many radio relics, including some of the best known 
and studied (e .g., the Coma cluster relic 1253+275 
and A bell 3667: IGiovannini et al.lll99ll ; iRottgering et al.1 
[1993), are located far from the cluster center, on the 
order of Mpcs. However, there exists a subclass of ra- 
dio relics which are located signi ficantly closer to the 
cluste r center, ~50 to 350 kpc. IGiovannini fc Ferettil 
(2004) classify these as relic sources near the first ranked 
galaxy. These radio sources share many of the same 
properties of classical radio relics including s teep spec- 
tra, fi lamentary structure, and polarization (Sle e et al.l 
l2T)0l . Abell 133 is one of the class of radio relics near 
the first ranked galaxy. A Chandra observation revealed 
a tongue-like feature in the X-ray emission comprised of 
cool gas [kT w 1.3 keV) which extended from the cD 
galaxy to the radio relic (|Fuiita et al.ll2002T ). A number 
of possible origins for the tongue were explored includ- 
ing the suggestion that the cool gas in the tongue was 
uplifted from the cluster core by the buoyantly rising ra- 
dio relic. A followup observation with XM M-Newton is 
consi stent with the uplifted bubble scenario (|Fuiita et al.l 
120041 ). However the detection of a cold front southeast 
of the cluster core and a weak shock in the core sug- 
gests that an unequal mass merger is responsible for the 
tongue feature. In either model, the radio plasma ap- 
pears to have originated from the active galactic nucleus 
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TABLE 1 

Properties of Detected Point Sources in Abell 13 Field 



Name 




R.A. a 


Dec. a 


Count Rate (10~ J s _1 ) 


CXOU J001316.6- 


193201 


00:13:16.67 


-19:32:01.2 


6.3±0.4 


CXOU J001319.0- 


193018 


00:13:19.09 


-19:30:18.5 


0.43±0.12 


CXOU J001320.0- 


192748 


00:13:20.06 


-19:27:48.6 


2.1±0.2 


CXOU J001321.9- 


192750 


00:13:21.98 


-19:27:50.7 


2.4±0.2 


CXOU J001326.0- 


193111 


00:13:26.02 


-19:31:11.4 


3.9±0.3 


CXOU J001335.7- 


192805 


00:13:35.74 


-19:28:05.5 


0.61±0.14 


CXOU J001339.9- 


192935 


00:13:39.94 


-19:29:35.8 


4.2±0.3 


CXOU J001341.1- 


192809 


00:13:41.15 


-19:28:09.5 


1.8±0.2 


CXOU J001342.5- 


193209 


00:13:42.58 


-19:32:09.0 


0.54±0.13 


CXOU J001344.6- 


192914 


00:13:44.60 


-19:29:14.1 


0.38±0.11 


CXOU J001344.7- 


193051 


00:13:44.70 


-19:30:51.0 


4.6±0.3 



a Units of right ascension are hours, minutes, and seconds, and units of declina- 
tion are degrees, arcminutes, and arcseconds. The 90% uncertainty in absolute 
positional accuracy of Chandra is CK'6. 



(AGN) in the central cD, rather than as a direct result 
of merger shock acceleration. 

In this paper, we present X-ray data obtained with 
the Chandra X-ray Observatory for the cluster Abell 
13, another unusual system with a radio relic near the 
first ranked galaxy. Abell 13 is an X-ray luminous 
(8 x 10 43 er g s" 1 in the .5-10 keV range) cluster at 
z = 0.0943 (jMazure et aLlll996l ). The relic has an ex- 
tremely steep (a w 4) radio spectrum, a linear po- 
larization of sal 2% at 1425 MHz, and spans up to 2' 
(jSlee et al.ll2001l ). No identification of a host galaxy as- 
sociated with the relic was possible, although the two 
brightest cluster members are candidates. Interestingly, 
these two galaxi es are separated by 2621 km s _1 in reces- 
sional velocity (|Quintana fc Ramirezl |1995[ ) . Similarly, 
from rcdshift dat a of the ESQ Nea rby Abell Clusters 
Survey (ENACS). lFadda et all (|1996t ) concluded that the 
radial velocity distribution of Abell 13 may also be in- 
terpreted as a bimodal one, with peaks at z = 0.0919 
(16 redshifts, vel.disp. 36lt^ km s" 1 ) and 0.0972 (21 
redshifts vel.disp. hlht™ km s" 1 ), or < cz >= 27570 
and 29160 km s _1 . We assume H Q = 70 km s" 1 Mpc -1 , 
= 0.3, and Q\ = 0.7 unless otherwise mentioned. At 
the distance of Abell 13, 1" corresponds to 1.75 kpc. 

2. SPATIAL ANALYSIS 

We observed Abell 13 on 2004 August 25 for 55 ks with 
Chandra us ing the Advanced C CD Imaging Spectrome- 
ter (ACIS: iGarmire et alJ 12003). The observation was 
performed so that both the cluster center and the radio 
relic were located on the back-illuminated S3 CCD. The 
"level 1" event file was processed using the CIAO v3.2 
data analysis package 10 following the standard reduction 
procedure. We removed the afterglow correction and cre- 
ated a new observation specific bad pixel file. The level 
1 event file was reprocessed to include the new bad pixel 
file and utilize the improved background event detection 
in VFAINT mode. We then filtered the data to remove 
hot pixels, bad columns, and events with grades 1, 3, 5, 
and 7. Data from the back-illuminated SI CCD were 
used to determine time intervals with high background 
using the tool lc_clean. We filtered the event file on 
good times, yielding a final exposure time of 50.7 ks. 
Background data were extracted from the blank sky ob- 

10 http://asc.harvard.edu/ciao/ 



servations compiled by M. Markevitch 11 . Our analysis is 
restricted to data from the S3 CCD. 

In Figure [TJ we show both the raw Chandra X-ray im- 
age of the 0.3-10 keV energy range, and an adaptively 
smoothed image of the X-ray data produced with the 
CIAO tool csmooth. The smoothed image has a mini- 
mum signal-to-noise ratio of 3 per smoothing beam, and 
has been corrected for exposure and background. The X- 
ray emission from Abell 13 shows complex structure in 
the core region. A knot of X-ray emission is also clearly 
detected 86" (150 kpc) to the northwest of the cluster's 
X-ray peak, along with a more diffuse extension of the 
emission in the same direction. 

We identify 11 point sources in our full band image 
(0.3-10.0 keV; see Table [J). This is consistent with 
the expected number of background sources from the 
Chan dra Deep Fields (jBrandt et al.ll2001_ iGiacconi et al.1 
2001). No point source was found at the center of the X- 
ray emission. Since the cluster emission will reduce our 
point source detection sensitivity, we also searched for 
point sources in just the hard band (4.0-10.0 keV) data. 
In the hard X-ray band, the diffuse emission from the 
cluster is weaker while the emission from AGN is stronger 
than at softer bands. Therefore, any AGNs should be 
more easily detected. We found no point sources within 
Q'.5 of the cluster center, defined to be the position of 
galaxy H, the brightest cluster galaxy. 

In Figure [2j we show the optical Digital Sky Survey 
(DSS) image 12 of the cluster with the X-ray contours 
overlaid. The X-ray contours are based on the image 
in Figure [1] right panel. We have l abeled the brightest 
galaxies following |Slee_et_aL| ((2001). The X -ray emis- 
sion is centered near the bright elliptical H (also 2MASX 
J00133853-1930007), the brightest cluster galaxy. The 
knot of emission to the NW is coincident with the ellip- 
tical galaxy F (also 2MASX J00133401-1929017), the 
second brightest cluster galaxy. The steepness of the X- 
ray gradient to the north and east is an artifact of the 
CCD edges. 

We are also interested in comparing the X-ray data 
with the radio emission from Abell 13. Figure [3] shows 
the X-ray data in red and the radio data in green. We 
find no enhancement of the X-ray emission near the ra- 
dio relic, which is different from what was seen for Abell 

11 See http://cxc.harvard.edu/cal/Acis/Cal_prods/bkgrnd/acisbg/ 

12 Available at http://archive.stsci.edu/dss/ 
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Fig. 1. — Left Panel: Raw Chandra X-ray image (0.3—10 keV) of the 5' X 4' region of Abell 13 encompassing the cluster center and radio 
relic positions. The image is uncorrected for exposure or background. North is up, and east is left. The image pixels are 0'.'492. Right 
Panel: Adaptively smoothed Chandra X-ray image of the same region, corrected for background, exposure, and vignetting. Notice the 
structure in the X-ray core as well as the knot of emission to the northwest of the core. In addition, the X-ray emission is extended to the 
northwest and west (towards the radio relic). The number of point sources is consistent with that expected from background sources. 

might expect the radio relic to extend if it were projected 
back to galaxy H. 

2.1. X-ray Emission from Galaxy F 

Enhanced X-ray emission has been s een in the large 
ellipt ical galaxies in the Coma cluster (Vikhlim n et al.l 
l2001f ). We compared the enhanced emission of galaxy 
F with the results from Coma. First, we determined 
the radial profile of the enhancement to determine its 
size and the appropriate source region for galaxy F. The 
emission was clearly more extended than the Chandra 
point spread function. We found that the enhancement 
measured «l'.'5-2" i n radius, or rou g hly 3 kpc at the 
distance of Abell 13. iVikhlinin et~all (|2001[) found that 
the compact cores associated with the brightest galaxies 
in the Coma cluster had a physical size of 3 kpc in radius. 

We attempted to extract a spectrum of the enhanced 
emission associated with galaxy F from a circular region 
with radius of 2". To account for both the cosmic X- 
ray background and the cluster emission, we used a lo- 
cal background extracted from an annulus centered on 
galaxy F with radii of 3" and 5". There are only 31 
counts in a 2" region around galaxy F, and we expect 
17.5 of these to be background. This is too few counts 
for spectral fitting or even to determine a useful color for 
the source. Even a basic color comparison cannot yield 
any constraint on the temperature of the galaxy gas. 

We used PIMMS to determine the flux from galaxy 
F. We measured a 0.5-8.0 keV count rate of 2.1 x 10~ 4 
cts s _1 and assumed a gas temperature of 1 keV for a 
Raymond-Smith plasma model. The largest elliptical 
galaxies have X-ray g as temperatures of sal keV (e.g., 
lO'Sullivan eFaT][2003l) . From PIMMS we find a 0.5-2.0 
keV unabsorbed flux of 8.2 x 10 -16 erg cm -2 s _1 , or a 



133 (jFuiita et alj 12002) . Given the low surface bright- 
ness of the X-ray emission, it is difficult to determine 
if there is a deficit of X-ray emission at the position of 
the radio relic, commonly referred to as an X-ray cavity 
or radio bubble. The absolute astrometric accuracy of 
the radio and X-ray maps is roughly 0'.'5 for each. The 
radio point source near the cluster center is not associ- 
ated with galaxy H. Instead, it seems coincident with a 
barely resolved pair of galaxies in the DSS and Two Mi- 
cron All Sky Survey (2MASS) images, listed in NED 13 as 
2MASX J00133712- 1930067. No velocity information is 
available for 2MASX J00133712- 1930067, making asso- 
ciation with the cluster uncertain. 

To enhance the structures in the cluster core, we fit- 
ted the observed cluster emission with a smooth ellip- 
tical isophotal model, using the IRAF/STSDAS task 
ellipse. We allowed the cllipticity, position angle, and 
intensities of each isophote to vary while keeping the cen- 
troids fixed to the position of galaxy H. The model was 
then subtracted from the data after multiplication of the 
model counts by 0.5 to avoid oversubtraction. In Fig- 
ure |H we show the subtracted image of the cluster in 
the region around the radio relic with the radio contours 
overlaid. No excess emission is associated with the radio 
relic. The cluster emission is faint near the radio relic, 
and continues to weaken at projected cluster-centric dis- 
tances greater than that of the relic in the same direction. 

More interesting structure is seen in the cluster center. 
Figure [5] shows the same image as Figure H] with a dif- 
ferent color scale to enhance the variations in the inner 
region of the cluster. The knot of emission in the up- 
per right is the second-brightest elliptical galaxy, F. The 
bright yellow knot of emission is coincident with galaxy 
H. There is a more diffuse bright extension of emission to 
the west of galaxy H. In addition, there is other structure 
in the cluster center perpendicular to the bright exten- 
sion. Interestingly, the emission extension to the west 
of galaxy H seems to point in the same direction as one 
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luminosity of 2 x 10 erg s . The galaxies in Coma had 
similar, although slightly larger (7.6-9.1 x 10 40 erg s _1 ) 
luminosities. We also calculated the mass of gas neces- 
sary to produce the excess emission from galaxy F. To 
do this, we first calculated the model normalization of a 
1 keV APEC plasma model that would produce the flux 
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Fig. 2. — X-ray brightness contours (0.3—10 keV band, logarithmically spaced by a factor of v2~) overlaid on the DSS opt ical image 
scann ed from a SERC-J survey plate taken in 1977 at the UK Schmidt telescope, Australia. The bright galaxies H and F from Sice ct al. 
II2001I) are labeled. The X-ray emission is centered near galaxy H. A knot of X-ray emission to the northwest of the cluster X-ray peak is 
consistent with the position of galaxy F. The steepness of the X-ray contours on the northern and eastern edges are artifacts of the CCD 
edge. 

3.1. Total Spectrum 

We first determined the integrated spectrum of the 
cluster using an elliptical region with semi-major and 
semi-minor axes of 2'5 and 1.'9, respectively, and position 
angle of 121° measured from north to east (see Figure^. 
This region was the largest elliptical region whose shape 
agreed with the outer cluster X-ray isophotes and fit en- 
tirely on the S3 CCD. We fit the data with an absorbed 
plasma model, using the APEC model in XSPEC. We fixed 
the redshift of the cluster to the mean re dshift of Abell 
13 fo und in the ENACS study [z = 0.0943: iMazure et al.l 
1996) and the equivalent hydrogen column density to 
the average Galactic val ue in the direction of Ab ell 13, 
N H = 2.0 x 10 20 cm~ 2 (jDickev fc Lockmadfl99(l . The 
best-fit model had a temperature kT — 6.0 ±0.3 keV and 
a metal abundance of 0.46 ± 0.12 solar (see Tableland 
Figure [7]) . 

We also performed a spectral study to determine if a 
cooling flow is present in the center of Abell 13, i.e. the 
peak in the X-ray emission which is also coincident with 
galaxy H. We extracted spectra in circular annuli out 
to 100" from the cluster center. Each spectrum was re- 
quired to have at least 1000 counts. The width of the 
annuli ranged from 15" in the cluster center to 5" at the 
periphery of the cluster. We fitted each spectrum with 
the same APEC model used to describe the total cluster 
emission, allowing only the temperature and normaliza- 
tion to vary (the metallicity was fixed at the best fit value 



measured from galaxy F. We then calculated the aver- 
age electron number density in the gas from the model 
normalization and assuming a spherical gas distribution 
with 2" radius. Translating the electron number den- 
sity into a gas mass, we find that the required gas mass 
is 4.2 x 10 7 Mq. This is comparable to the gas masses 
found in the Coma galaxies and other galaxies in nearby 
clusters ()Vikhlinin et al.ll200H ISun et al.ll2007f) . 

Given the similarity in size and luminosity, we may ex- 
pect that the cause of the X-ray enhancement in galaxy 
F is the same as th at found for the Coma galaxies 
(jVikhlinin et al.|[200l . The gas in the Coma galaxies 
is of low temperature (1-2 k eV) and is likely t he rem ains 
of the galactic X-ray halos. Vikhlin m et alj (|2001f ) sug- 
gest that the cool gas is in pressure equilibrium with the 
hotter, intercluster medium. 

3. SPECTRAL ANALYSIS 

Spectra were extracted from the level 2 event file, 
excluding point sources within the regions of interest. 
Background spectra were extracted from the blank-sky 
background file. For each spectrum, weighted response 
files (RMFs and ARFs) were generated using the CIAO 
tools mkacisrmf and mkwarf , respectively. Finally, spec- 
tra were grouped to obtain a minimum of 25 counts per 
bin. Spectral fitting was primarily performed in XSPEC 
vll.0, except for the temperature map spectral fits (see 
§EJ below). 
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Fig. 3. — Radio image at 1.4 GHz (Slcc ct al. 2001) overlaid in green the adaptively smoothed X-ray image of Abell 13 in red. North is 
up, and east is le ft. The region covered is 5'.l X 3.'8. The radio image shows only emission detected at greater than 2.5cr as determined by 
Slcc ct al. (2001). The three compact red spots are X-ray point sources unrelated to the cluster. 

est to the count-weighted average position of the events 
in each region. The spectra were grouped to have 25 
counts per bin and were fitted in ISIS using the APEC 
thermal model. The Galactic hydrogen column density, 
redshift, and abundance were fixed to the values found 
for the full cluster spectrum (§ 13.11 above). Only the 
temperature and normalization of the APEC model were 
allowed to vary. In general, the errors on the tempera- 
ture map are about 20% in the temperature. Although, 
for pixels with higher best-fit temperatures or those near 
the edge of the map, the errors are larger, in the range 
of 50-70%. This is because of the low high-energy sen- 
sitivity of Chandra and the lower statistics of the outer 
regions. 

The temperature map shown in Figure[8]is colored such 
that the coolest regions (kT « 4 keV) are red, and the 
hottest regions (kT » 12 keV) purple. The intensity 
contours for the X-ray and radio emission from Abell 13 
are overlaid. We note that small differences in tempera- 
ture from pixel to pixel are not statistically meaningful 
given the errors in the fits (an average error of 1.2 keV 
at 6.0 keV). Since the fitted region size is larger than a 
single map pixel, there is overlap of the fitted regions. 

Of particular interest are the lack of very cool gas in 
the cluster center and the presence of cool gas associated 
with the radio relic region. The temperature map shows 
a large concentration of lower temperature gas associated 
with the radio relic. Below f H3.3p we perform a detailed 
spectral study of this gas to better constrain its prop- 
erties. There may also be some hotter gas (purple/blue 



found above). No statistically significant variations in 
temperature were found. This result is consistent with 
the temperature map determined in § 13.21 

3.2. Temperature Map 

In order to determine the spectral distribution 
in the cluster emission, we created an adaptively 
binned temperature map of Abell 13, shown in Fig- 
ure [8] The temperature map was produced using ISIS 
(jHouck &: Deni cola 2000]). A spectrum at each map pixel 
(1 map pixel = 16.3 ACIS pixels = 8'.'0) was extracted 
from the event data (excluding point sources) in a square 
region centered on the pixel. The region size was deter- 
mined by requiring at least 800 net (i.e. background 
subtracted) counts in the 0.4-7.0 keV energy range, but 
with a maximum box size of 183x183 ACIS pixels, or 
90"x90". The maximum box size was only reached in 
the outskirts of the cluster. The total counts in each 
region were between 1.5 and 2 times the 800 count min- 
imum to account for background subtraction. For each 
spectrum, a background spectrum was extracted from 
the blank sky background event file in the same region. 
RMF and ARF files were constructed for specific loca- 
tions on the CCD, with a spacing of 32 and 50 pixels, 
respectively. This was done to account for variations 
in the response across the CCD, without having to cre- 
ate individual RMF and ARF files for each region which 
would be computationally demanding. For a given re- 
gion, the RMF and ARF used to fit the spectrum were 
chosen by selecting the response files with positions clos- 
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Fig. 4. — Residual X-ray emission from the cluster (color scale from blue=weak to white=strong) after partial subtraction of a smooth 
elliptical isophotal model for the entire cluster. The region shown here encompasses the radio relic and has a size of 4'6 X 2.'8. The color 
scale w as chosen to emphasize the outer regions of the cluster near the radio relic. Overlayed are the radio contours as given in lSlee et al.l 
(|2001f ) with contour levels of 43, 86, 173, 303 and 389 /^Jy beam -1 , which correspond to 2.5, 5, 10, 18 and 23 a. The two brightest galaxies 
H and F are indicated by the Xs. No excess or clearly associated deficit in X-ray emission is found coincident with the radio relic. Along 
the eastern tail of the radio emission there appears to be a sharp drop in the residual X-ray emission. However, this drop is not very strong 
in the total X-ray surface brightness. There is a strong gradient in the southern part of the cluster, but at a larger radius and further south. 

region) associated with the steep gradient in the X-ray 
surface brightness contours on the southern edge of the 
cluster (see e.g., the X-ray contours in Figure [5]) but the 
temperature errors are large in this region. To try and 
confirm this result, we extracted spectra on the northern 
and southern sides of the X-ray gradient. No signifi- 
cant difference in the best-fit temperatures was found, 
although both regions had large errors. Higher count 
spectra are required to verify the increased gas temper- 
ature along the cluster's southern edge. 

3.3. X-ray Spectral Study of the Radio Relic Region 

From the temperature map, we see that the X-ray 
emission associated with the radio relic is colder than the 
average cluster temperature. This immediately rules out 
a shock origin as we would expect the emission to be hot- 
ter in that case. To determine the nature of this cooler 
emission, we performed a more rigorous X-ray spectral 
study of the radio relic region. We fit the spectra from 
a 2.'0xlf4 (major/minor axis) elliptical region surround- 
ing the radio relic (Figure [5]) . This region was chosen 
to encompass the outer portion of the radio relic that is 
associated with the cooler gas in the temperature map. 
This region has the advantage that it is a region of lower 
surface brightness for the hot gas component. If we add 
in the radio relic region to the east (the tail leading back 
to galaxy H), the spectrum has a much larger hot gas 
component and is less sensitive to a cool gas component. 

The data were fit with different spectral models, 
including a single temperature gas model, a two- 
temperature gas model, and a gas + powerlaw emission 
model. The results are given in Table [3] In all fits we 
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Fig. 5. — Residual emission from center of Abell 13 after partial 
subtraction of a smooth elliptical isophotal model. The figure spans 
a region of 2'4 X 2.'2. The color scale was chosen to emphasize the 
inner region of the cluster near galaxy H. North is up and east is 
left. The bright region marked with an X near the center coincides 
with galaxy H. Notice the extension of emission to the west (in the 
direction of the radio relic) and the excess emission perpendicular 
to that direction. The X at the top right corner marks the position 
of galaxy F. 
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Fig. 6.— Raw Chandra X-ray image (0.3-10 keV) of Abell 13 
encompassing the cluste r center and radio relic positions. Overlaid 
are the radio contours of Slec ct al. (2001). The straight lines mark 
the coverage of the VLA data compared to the Chandra data. The 
large ellipse marks the region of the data use d to determine the 
spectral model of the full cluster emission (§ 13.11 1. The smaller 
ellipse on the right marks the region of the data used to measure 
the spectrum of the X-ray gas coincident with the radio relic (§ 13.31 1. 
The region covered is 6.'6 X 4.'7. 



TABLE 2 

Best-Fit APEC Model Parameters for 
Cluster 3 - 



Parameter Value 



AT H (cm~ 2 ) 2.0 x 10 20b 

kT (keV) 6.0±0.3 
Abundance (solar) 0.46±0.12 

Norm (2.98±0.08)xl0" 3 
X 2 /dof 316.45/297 



a All errors quoted at the 90%-confidence 
level. 

b Column density fixed at Galactic value. 

have frozen the metallicity of the thermal components to 
that of the general cluster emission. If allowed to vary, 
the metallicity is unconstrained. The single temperature 
gas model, while a decent fit to the data, does not repro- 
duce the excess emission in the data near 0.8 keV (see 
Figure O. The two-temperature gas model is the best fit 
to the data. One of the components is consistent with 
the overall cluster temperature, while the other compo- 
nent has a temperature of «0.9 keV and accounts for 
the soft excess seen in the data. The gas + powerlaw 
model is a worse fit, with the powerlaw component re- 
placing the hotter gas component. The two-temperature 
gas model is the most appropriate, as it would include 
both a lower temperature gas associated with the relic 
as well as hotter cluster emission seen in projection. 

4. DISCUSSION 

We have presented the Chandra data of the galaxy clus- 
ter Abell 13. The X- ray emission is centered on galaxy 
H of ISlee et all (|2001h . the brightest cluster galaxy. The 
spatial structure near the core is complicated, with a 
bright tail of excess emission to the west of galaxy H, 
and a fainter excess to the north. No point source is 
associated with galaxy H and no cooling flow is found. 



1.000F 
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Fig. 7. — Upper panel: Count spectrum of the X-ray emission 
from Abell 13 and the best fit model. Lower panel: Ratio of the 
data and best-fit model (data/model). The spectrum was extracted 
from a 2.'5 X 1.'9 elliptical region encompassing the largest elliptical 
isophote contained fully on the S3 CCD. The best-fit model has a 
temperature kT = 6.0 ± 0.3 keV and a metal abundance of 0.46 ± 
0.12 solar. 

We also found excess X-ray emission associated with the 
second-brightest elliptical galaxy, F, in the cluster. The 
mass of gas associated with galaxy F is consistent with 
other bright elliptica l galaxies in clu ster environments 
(jVikhlinin et al.ll200HlSun et al.ll2007| ). In general, there 
is an extension of the X-ray emission to the northwest 
with respect to the X-ray peak coincident with galaxy 
H. Finally, no excess or significant deficit in emission is 
found coincident with the unusual radio relic in Abell 13, 
although the spectral properties of the gas in this region 
are different from those in other areas of the cluster. We 
found a soft excess in the spectrum of gas coincident with 
the radio relic, which can be well fit by a «0.9 keV ther- 
mal model, suggesting the presence of cooler gas there. 

The redshift distribution of the galaxies in Abell 13 
can b e described by a bimodal distribution ([Fadda et al.l 
1996). The subclumps are not distinguishable in their lo- 
cation on the sky. Interestingly, the two brightest galax- 
ies, H and F, are associated with different peaks in the 
redshift distribution and rather than being associated 
with the peak velocity, either of the cluster as a whole or 
of either subclump, the velocities of galaxies H and F are 
on the very edges of the velocity distribution. Galaxy F 
is on the lower edge for the low-velocity subcluster and 
galaxy H on the higher edge of the high- velocity subclus- 
ter. This, combined with the complicated structure of 
the X-ray emission, points to an ongoing merger within 
the cluster and that galaxies H and F have a significant 
velocity relative to the general cluster distribution. The 
spatial distribution of the galaxies in Abell 13 shows a 
northwest-southeast extension, paralleling the X-ray dis- 
tribution, but extending to larger scales. 

In the following discussion, we assume that the X-ray 
gas associated with the radio relic has two components, 
a hot (6 keV) gas component due to the general cluster 
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Fig. 8. — Temperature map of Abell 13. Red signifies regions of lower temperature, purple higher temperature. Overlaid are the X-ray 
and radio surface brightness contours. No cooling flow is detected. The coolest gas is to the west of the peak in the X-ray emission, and 
apparently associated with the radio relic. 



emission, and a cooler (0.9 keV) gas component local to 
the relic. This model is based on our best-fit spectrum 
found in § 13.31 We assume that the two components are 
in press ure equilib r ium. 

From IsTee et al.l (|2001h . the minimum pressure of the 
radio relic in Abell 13 is 1.4 x 10~ 12 dyne cm -2 . Their 
estimate of the thermal pressure of the hot X-ray gas was 
very rough and based solely on scaling the properties of 
Abell 85 since no imaging observation of Abell 13 was 
available. Our Chandra data allow for a more accurate 
estimate of the thermal pressure at the location of the 
relic. However, due to an unfortunate placement of the 
cluster on the CCD, the only regions which extended out 
to a distance larger than the relic were in the direction 
of the relic and to the south. Since the southern edge 
of the cluster shows a steep gradient in the X-ray emis- 
sion, we did not consider this region in determining the 
surface brightness distribution. Instead, we determined 
the surface brightness of the cluster in a conical region 
extending from the cluster center out to 4fl and encom- 
passing position angles of 240-312° measured from the 
north to the east. This region includes the radio relic and 
associated cool gas which will affect our results; however, 
while not a perfect measure of the properties of the X- 



ray gas, it r epresents a signi ficant improvement on the 
estimates of ISlee et al.1 (l200l . A larger field of view ob- 
servation would be better suited to address the surface 
brightness and electron density distributions of the full 
cluster emission. 

Using the surface brightness distribution, we esti- 
mate the electron density distribution of the hot com- 
ponent foll o wing the spherical deprojection method of 
iKriss et al.l (fl9 83 ) . We assumed a single temperature of 
the hot gas of 6.0 keV and used the APEC thermal model 
to convert from counts to electron density. From this, 
we find an electron density at the relic of 1 x 10~ 3 cm" 3 . 
Given this and the average cluster temperature (6.0 keV) 
determined in § 13. 1( we find that the thermal pressure of 
the surrounding hot gas at the location of the radio relic 
is 1.6x 10 -11 dyne cm~ 2 . This value is larger than the es- 
timate for the pressure of the radio plasma, as has been 
seen in other X-ray cluster radio bubbles (|Dunn et al.1 
2005). We note that if the relic were at a larger radius 
than the projected radius, the electron density and ther- 
mal pressure would be lower, bringing it closer to the 
radio plasma pressure estimate. In addition, in the case 
of Abell 13, if the values assumed bv lSlee et alj (|2001[ ) for 
(1) the ratio of masses in heavy particles and electrons; 
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TABLE 3 

Best-Fit X-ray Spectral Model Parameters for Radio RELic a 
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a All errors quoted at the 90%-confidence level. All models were fit including Galactic 
absorption of 2.0 X 10 20 cm -2 . Abundances for gas emission models were fixed to the best- 
fit value from the cluster spectrum, 0.46 of solar. The model parameters Normi and Norm2 
are the values of the first and second model component normalizations. kT 2 or T is the value 
of the second model component temperature if APEC or photon index if powerlaw. 
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Fig. 9. — Count spectra of the X-ray emission associated with 
the radio relic, fit with different models. We have fit the spectrum 
with a single temperature APEC model (top panel), two APEC mod- 
els (middle panel), and a powerlaw + APEC model (bottom panel). 
Notice the soft excess at 0.8 keV in the single temperature model 
(top) but not the other spectra. The peak at ~1.1 keV is too nar- 
row to be astrophysical and is likely strengthened by binning. At 
other binnings, it is not as significant. 

(2) the plasma filling factor; (3) the low frequency spec- 
tral index, and (4) the angle between the line of sight 
and field direction are given more realistic values then 
we may increase the radio plasma pressure to that of the 
confining X-ray gas. 

4.1. Origin of the Radio Relic 

It is unlikely that the radio relic in Abell 13 is a periph- 
eral cluster relic associated with a cluster merger shock. 



First, there is no evidence in the Chandra image or spec- 
tra for any shock feature near the relic. Second, the radio 
image shows a long filament which extends almost all the 
way back to the brightest elliptical galaxy H (Figure [3]). 
Thus, we assume that the radio emitting plasma origi- 
nated in an AGN at the center of galaxy H. B ased on 
the classification system of lKempner et all (|2004l ) , which 
looks at the physical properties and proposed origins of 
extended cluster radio sources, we would classify Abell 
13 as an AGN relic system. 

We propose two explanations for the multi- wavelength 
properties of Abell 13 . The first i s the b uoyant bubble 
scenario sugg ested bv lFuiita etafl ((2002) for Abell 133. 
The radio source may have uplifted the cool gas from 
the core of Abell 13. In this scenario, galaxy H would 
have been an active radio emitter in the past. At some 
later time, the radio plasma would have risen by buoy- 
ancy, and entrained the cool gas in the center, uplifting it 
to a larger distance from the cluster center. This might 
also remove the source of fuel for the central AGN. This 
would explain the lack of cool gas in the cluster center 
and the fact that currently , galaxy H is n either a ra- 
dio (< 0.1 mJy at 1.4 GHz; Isiee et aLll2001h nor a hard 
X-ray source. To estimate the time scale for the cool 
gas and radio source to rise by buoyancy from the posi- 
tion of galaxy H, we assume that the cool gas and radio 
source move at a velocity roughly 1/3 the sound speed, 
c s , of the hot gas component. For a 6.0 keV gas, c s 
is 980 km s _1 . The projected separation between the 
relic and galaxy H is 2'2-3.'l. This yields a buoyancy 
time scale of 7-10 x 10 s yr, roughly consiste nt with the 
relic age determined through spectral fitting (|Slee et al.1 
|2001| ). To estimate the mass of cool gas uplifted, we take 
the volume to be that of a prolate spheroid with ma- 
jor and minor axes equal to that of the region used to 
extract the spectrum of the radio relic. From the two- 
temperature spectral model, we find that the enclosed 
mass of cool gas is ~ 2 x 10 10 M Q . This is a factor of 10 
greater than the mass of gas c ontained in the to ngue of 
X-ray gas found in Abell 133 (jFuiita et al.ll2002f) . How- 
ever, we note that the electron density for the cool gas 
is ~ 5 x 10~ 4 cm~ 3 , lower than the hotter cluster gas 
(particularly the gas towards the center of the cluster) 
making buoyant uplift a possibility. A cool gas mass of 
~ 2 x 10 10 M Q implies that previously Abell 13 would 
have had a mass deposition rate of ~ 2 M Q yr -1 assum- 
ing a cooling time of 10 10 yr. This is on the low end of 
the distributio n of known coolin g flow mass deposition 
rates (see e.g., IWhite et al.lll997l ). This is only an esti- 
mate of the mass of cool gas to determine if the value 
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is realistic. The temperature map shows additional cool 
regions to the north and east of the radio relic which if 
taken into account would increase the mass of cool gas. 
However, the radio relic is filamentary in structure and 
it is possible that the cool gas also resides in filaments 
rather than a single big blob, which would decrease the 
gas mass estimate. Since it is difficult to determine the 
extent of either of these situations, we merely suggest 
that the estimated amount of cool gas in Abell 13 seems 
reasonable for a galaxy cluster. 

The lower temperature of the gas near the radio source 
might be due both to the lower entropy of the gas in the 
cool core of the cluster, plus the effects of adiabatic ex- 
pansion as this gas was moved out to a lower-pressure 
region. From our above estimate of the electron density 
distribution of the cluster gas, we find that the electron 
density at the cluster center is «6 times that at the radio 
relic. The pressure difference between the two locations 
should also be 6 times assuming a constant hot gas tem- 
perature. We assume that the cool gas component is 
always in pressure equilibrium with the hot gas compo- 
nent that dominates the cluster emission. As a result of 
adiabatic expansion, the temperature of the cool gas in 
its current location, T 2 is related to its initial tempera- 
ture, Ti, by 

where Pi is the pressure in the cluster core and Pi is 
the pressure at the location of the radio relic, and the 
adiabatic index is 7 = 5/3. Observationally, we found 
Pi/ Pi ~ 6, giving T 2 « 0.5 7\ for adiabatic expansion. 
This difference alone cannot account for the cooling of 
the gas from the average cluster temperature of «6 keV 
to its current temperature of w 0.9 keV. Therefore, the 
cool gas associated with the radio relic must have been 
cooler than the average cluster temperature to begin 
with. We can estimate the initial cool gas temperature 
assuming that adiabatic expansion was the only heat- 
ing or cooling mechanism as the cool gas was uplifted. 
This requires that the initial temperature of the cool gas 
was about T\ w 2 keV when this gas was located at the 
cluster core. This is a reasonable temperature for a cool- 
ing core. If the relic were at a larger radius than the 
projected radius, then we would expect the electron den- 
sity to be lower and therefore the adiabatic cooling to be 
stronger. This would make the initial temperature of the 
cool gas higher, possibly more in line with the average 
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cluster temperature. 

Our second proposed explanation is that a dynamical 
event is the origin of the unusual radio relic and disturbed 
X-ray morphology. Both the X-ray and optical proper- 
ties (radial velocity distribution) suggest that Abell 13 
is undergoing a merger event. The X-ray extension from 
galaxy H (to the west) points back to the radio relic. This 
may imply that galaxy H is moving to the east and that 
the radio relic is at the former position of the galaxy. 
(We note that the extension may also favor the buoy- 
ant uplift scenario.) We propose that before the merger 
event, galaxy H had a cool core and showed AGN ac- 
tivity. As the merger progressed, the cold core and ra- 
dio lobes were displaced (or stripped) from the galaxy 
as a result of ram pressure from the gas in the other 
subcluster, whic h did not a f fect th e stars and dark mat- 
ter in galaxy H. ISlee et ail (|2001h found that the travel 
time from galaxy H was comparable to the estimated 
relic age. The long filament of radio emission leading 
back almost to galaxy H (Figure [3]) suggests that the 
AGN remained active during the initial phases of this 
stripping process. The concentration of X-ray emission 
around galaxy H (Figure suggests that the stripping 
process is not complete; there is still gas around galaxy 
H, but it is relatively hot, perhaps due to heating pro- 
cesses associated with the merger. The Chandra data do 
not conclusively determine which of these two models is 
correct, or strongly reject either. 
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